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Abstract. A review is presented of the fireball model of 7-ray bursts (GRBs), and 
of the production in GRB fireballs of high energy protons and neutrinos. Constraints 
imposed on the model by recent afterglow observations, which support the association 
of GRB and ultra-high energy cosmic-ray (UHECR) sources, are discussed. Predictions 
of the GRB model for UHECR production, which can be tested with planned large area 
UHECR detectors and with planned high energy neutrino telescopes, are reviewed. 
PACS numbers: 98.70.Rz, 98.70.Sa, 96.40. Tv, 14.60. Pq 



1 Introduction and summary 

The widely accepted interpretation of the phenomenology of 7-ray bursts (GRBs), 
bursts of 0.1 MeV-1 MeV photons lasting for a few seconds (see for review), 
is that the observable effects are due to the dissipation of the kinetic energy 
of a relativistically expanding wind, a "fireball," whose primal cause is not yet 
known (see |6^,Q for reviews). The recent detection of "afterglows," delayed low 
energy (X-ray to radio) emission of GRBs (see for review), confirmed the 
cosmological origin of the bursts, through the rcdshift determination of several 
GRB host-galaxies, and confirmed standard model predictions of afterglows that 
result from the collision of an expanding fireball with its surrounding medium 
(see 1^ for review). In this review, the production in GRB fireballs of 7-rays, 
high-energy cosmic-rays and neutrinos is discussed in the light of recent GRB 
and ultra-high-energy cosmic-ray observations. 

The fireball model is described in detail in §|^. We do not discuss in this 
section the issue of GRB progenitors, i.e. the underlying sources producing the 
relativistic fireballs. At present, the two leading progenitor scenarios are collapses 



of massive stars [36,Sl|, and mergers of compact objects [42 79[. As explained in 
§1^, the evolution of the fireball and the emission of 7-rays and afterglow radiation 
(on time scale of a day and longer) are largely independent of the nature of the 
progenitor. Thus, although present observations provide stringent constraints 
on the fireball model, the underlying progenitors remain unknown (e.g. p^ ; see 
[ |6lU70| | for discussion) . In §||, constraints imposed on the fireball model by recent 
afterglow observations are discussed, which are of importance for high energy 
particle production. 
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The association of GRBs and ultra-high energy cosmic-rays (UHECRs) is 
discussed in Recent afterglow observations strengthen the evidence for GRB 
and UHECR association, which is based on two key points (see [108| for recent 
review). First, the constraints imposed on fireball model parameters by recent 
observations imply that acceleration of protons is possible to energy higher than 
previously assumed, ~ 10^^ eV. Second, the inferred local (z — 0) GRB energy 
generation rate of 7-rays, ^ lO^'^erg/Mpc'^yr, is remarkably similar to the local 
generation rate of UHECRs implied by cosmic-ray observations. 

The GRB model for UHECR production makes unique predictions, that may 
be tested with operating and planned large area UHECR detectors p3|,^5|jl0^j97|| . 
These predictions are described in f|5[ In particular, a critical energy is predicted 
to exist, 10^°eV < Ec < 4 x lO^^eV, above which a few sources produce most 
of the UHECR flux, and the observed spectra of these sources is predicted to 
be narrow, Z\e/e ^ 1: the bright sources at high energy should be absent in 
UHECRs of much lower energy, since particles take longer to arrive the lower 
their energy. If the sources predicted by this model are detected by planned 
large area UHECR detectors, this would not only confirm the GRB model for 
UHECR production, but will also provide constraints on the unknown structure 
and strength of the inter-galactic magnetic field. 

We note, that the AGASA experiment has recently reported the presence of 
one triplet and 3 doublets of UHECR events above 4 x lO^^eV, with angular 
separations (within each group) < 2.5°, roughly consistent with the measure- 
ment error [Q. The probability that these multiplets are chance coincidences 
(as opposed to being produced by point sources) is ~ 1%. Therefore, this ob- 
servation favors the bursting source model, although more data are needed to 
confirm it. Testing the predictions of the fireball model for UHECR production 
would require an exposure 10 times larger than that of present experiments. 
Such increase is expected to be provided by the HiRes a nd Auger ||25,103| 
detectors, and by the proposed Telescope Array detector ||97[. 

Predictions for the emission of high energy neutrinos from GRB fireballs are 
discussed in §|| Implications for planned high energy neutrino telescopes (the 
ICECUBE extension of AMANDA, ANTARES, NESTOR; see |6| for review) 
are discussed in detail in §6^. It is shown that the predicted flux of > lO^"* eV 
neutrinos may be detectable by Cerenkov neutrino telescopes while the flux 
above lO'^^ eV may be detectable by large air-shower detectors p9| , |65| , p^ . De- 
tection of the predicted neutrino signal will confirm the GRB fireball model for 
UHECR production and may allow to discriminate between different fireball 
progenitor scenarios. Moreover, a detection of even a handful of neutrino events 
correlated with GRBs will allow to test for neutrino properties, e.g. flavor oscil- 
lation and coupling to gravity, with accuracy many orders of magnitude better 
than currently possible. 



2 The Fireball model 
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2.1 Relativistic expansion 

General phenomenological considerations, based on 7-ray observations, indicate 
that, regardless of the nature of the underlying sources, GRBs are produced 
by the dissipation of the kinetic energy of a relativistic expanding fireball. The 
rapid rise time and short duration, ^ 1 ms, observed in some bursts p^ , ^ 
imply that the sources are compact, with a linear scale comparable to a light- 
ms, ro ~ 10'' cm. The high 7-ray luminosity implied by cosmological distances, 

~ lO^^erg s~^, then results in a very high optical depth to pair creation. The 
energy of observed 7-ray photons is above the threshold for pair production. 
The number density of photons at the source n~f is approximately given by 

= A-KrQcn^e, where e ~ IMcV is the characteristic photon energy. Using 
ro ^ lO^cm, the optical depth for pair production at the source is 

T^^ ~ ron^ar ~ .^ll^ ^ xO^^, (1) 
47rroce 

where is the Thomson cross section. 

The high optical depth implies that a thermal plasma of photons, electrons 
and positrons is created, a "fireball," which then expands and accelerates to 
relativistic velocities |7^,^ ■ The optical depth is reduced by relativistic expan- 
sion of the source: If the source expands with a Lorentz factor F, the energy of 
photons in the source frame is smaller by a factor F compared to that in the 
observer frame, and most photons may therefore be below the pair production 
threshold. 

A lower limit for F may be obtained in the following way [p7p|]. The GRB 
photon spectrum is well fitted in the Burst and Transient Source Experiment 
(BATSE) detectors range, 20 keV to 2 MeV by a combination of two power- 
laws, dn^/de^ oc e^^ with different values of [3 at low and high energy 
Here, dn^/de^ is the number of photons per unit photon energy. The break 
energy (where /3 changes) in the observer frame is typically t^b ~ IMeV, with 
/3 ~ 1 at energies below the break and /? ~ 2 above the break. In several cases, 
the spectrum was observed to extend to energies > 100 MeV (9|,^. Consider 
then a high energy test photon, with observed energy tt, trying to escape the 
relativistically expanding source. Assuming that in the source rest frame the 
photon distribution is isotropic, and that the spectrum of high energy photons 
follows dn^/de^ oc e~^, the mean free path for pair production (in the source 
rest frame) for a photon of energy = tt/l (i^i the source rest frame) is 



'77 



Here, eth(eti 0) is the minimum energy of photons that may produce pairs inter- 
acting with the test photon, given by ethet(l — cos 61) > 2(meC^)^ {9 is the angle 
between the photons' momentum vectors), t/^ is the photon energy density (in 
the range corresponding to the observed BATSE range) in the source rest-frame, 
given by L~f — Airr'^j'^cU-y. Note , that we have used a constant cross section. 
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3crT/16, above the threshold eth- The cross section drops as log(e)/e for e ^ eth', 
however, since the number density of photons drops rapidly with energy, this 
does not introduce a large correction to l-^-y. 

The source size constraint implied by the variability time is modified for a 
relativistically expanding source. Since in the observer frame almost all photons 
propagate at a direction making an angle < 1 /F with respect to the expansion 
direction, radiation seen by a distant observer originates from a conical section 
of the source around the source-observer line of sight, with opening angle ~ 1/^- 
Photons which are emitted from the edge of the cone are delayed, compared to 
those emitted on the line of sight, by r/2r^c. Thus, the constraint on source size 
implied by variability on time scale At is 

r - 2r^cAt. (3) 

The time r/c required for significant source expansion corresponds to comoving 
time (measured in the source frame) tco. ~ r/Fc. The two-photon collision rate at 
the source frame is t~]j = c/l^^. Thus, the source optical depth to pair production 
is Tj-y = tco./^77 ^ r/ri^~j. Using Eqs. (||) and (^) we have 

1 arL^et 



12877 c2(mec2)2rMr 

Requiring r^^ < 1 at et we obtain a lower limit for F, 

1/6 



F > 250 



L-y 52 ( — ) Arl] , (5) 

. ^' VlOOMeV/ ' 



where = lO^^^,^ ^gerg/s and At = 10 ^At-2 s. 



2.2 Fireball evolution 

As the fireball expands it cools, the photon temperature F^ in the fireball frame 
decreases, and most pairs annihilate. Once the pair density is sufficiently low, 
photons may escape. However, if the observed radiation is due to photons escap- 
ing the fireball as it becomes optically thin, two problems arise. First, the photon 
spectrum is quasi-thermal, in contrast with observations. Second, the source size, 
ro ~ 10^ cm, and the total energy emitted in 7-rays, ~ 10^'^ erg, suggest that 
the underlying energy source is related to the gravitational collapse of ^ IMq 
object. Thus, the plasma is expected to be "loaded" with baryons which may 
be injected with the radiation or present in the atmosphere surrounding the 
source. A small baryonic load, > 1O~^M0, increases the optical depth due to 
Thomson scattering on electrons associated with the "loading" protons, so that 
most of the radiation energy is converted to kinetic energy of the relativisti- 
cally expanding baryons before the plasma becomes optically thin [^^^. To 
overcome both problems it was proposed |Q that the observed burst is pro- 
duced once the kinetic energy of the ultra-relativistic ejecta is re-randomized by 
some dissipation process at large radius, beyond the Thomson photosphere, and 
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then radiated as 7-rays. Collision of the relativistic baryons with the inter-stellar 



medium and internal collisions within the ejecta itself [[78 83 7l|], were pro- 
posed as possible dissipation processes. Most GRBs show variability on time 
scales much shorter than (typically one hundredth of) the total GRB duration. 
Such variability is hard to explain in models where the energy dissipation is 
due to external shocks ]117] , p8[ . Thus, it is believed that internal collisions are 
responsible for the emission of 7-rays. 

Let us first consider the case where the energy release from the source is 
"instantaneous," i.e. on a time scale rg/c. We assume that most of the energy 
is released in the form of photons, i.e. that the fraction of energy carried by 
baryon rest mass M satisfies rj~^ = Mc^/E <g; 1. The initial thickness of the 
fireball shell is tq. Since the plasma accelerates to relativistic velocity, all fluid 
elements move with velocity close to c, and the shell thickness remains constant 
at rg (this breaks at very late time, as discussed below). We are interested in 
the stage where the optical depth (due to pairs and/or electrons associated with 
baryons) is high, but only a small fraction of the energy is carried by pairs. 

The entropy of a fluid component with zero chemical potential is S = V{e + 
p)lT, where e, p and V are the (rest frame) energy density, pressure and vol- 
ume. For the photons p = e/3 oc T^. Since initially both the rest mass and 
thermal energy of baryons is negligible, the entropy is provided by the photons. 
Conservation of entropy implies 

r^l{r)rQT^{r) ^ Constant, (6) 

and conservation of energy implies 

r'^-i{r)ra-i{r)T^{r) = Constant. (7) 

Here 7(7") is the shell Lorentz factor. Combining (^ and (|^) we find 

7(r) cx r, Tj{r) cx r^^, n oc r^^, (8) 

where n is the rest frame (comoving) baryon number density. 

As the shell accelerates the baryon kinetic energy, ^Mc^, increases. It be- 
comes comparable to the total fireball energy when 7 ^ 77, at radius rf = r]rQ. 
At this radius most of the energy of the fireball is carried by the baryon ki- 
netic energy, and the shell does not accelerate further. Eq. (|^) describing energy 
conservation is replaced with 7 = Constant. Eq. (|^), however, still holds. Eq. 
(^ may be written as T^/nT^ = Constant (constant entropy per baryon). This 
implies that the ratio of radiation energy density to thermal energy density as- 
sociated with the baryons is r independent. Thus, the thermal energy associated 
with the baryons may be neglected at all times, and Eq. (^ holds also for the 
stage where most of the fireball energy is carried by the baryon kinetic energy. 
Thus, for r > rf we have 

7(r) — r rj, T oc r^^^^, n cx . (9) 

Let us consider now the case of extended emission from the source, on time 
scale ^ tq/c. In this case, the source continuously emits energy at a rate L, and 
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the energy emission is accompanied by mass loss rate M = L/r](?. For r <rf the 
fluid energy density is relativistic, aT^/nm,pC^ = rjro/r, and the speed of sound 
is ~ c. The time it takes the shell at radius r to expand significantly is r/c in 
the observer frame, corresponding to tco. ~ r/'yc in the shell frame. During this 
time sound waves can travel a distance cr/^c in the shell frame, corresponding 
to = r/{r/ro)^ = {ro/r)ro in the observer frame. This implies that at the 
early stages of evolution, r ~ ro, sound waves had enough time to smooth out 
spatial fluctuations in the fireball over a scale ro, but that regions separated 
by Ar > ro can not interact with each other. Thus, if the emission extends 
over a time Tgrb tq/c, a fireball of thickness cTgrb ^ ''o would be formed, 
which would expand as a collection of independent, roughly uniform, sub-shells 
of thickness ro- Each sub-shell would reach a final Lorentz factor Ff, which may 
vary between sub-shells. This implies that different sub-shells may have velocities 
differing by Av c/2rf' , where 77 is some typical value representative of the entire 
fireball. Different shells emitted at times differing by At, tq/c < At < Tgrb, 
may therefore collide with each other after a time tc ~ cAt/Av, i.e. at a radius 

n « 2r'^cAt = 6 X 10^^rlr,At-2 cm, (10) 

where F = lO^'^/lj.s. The minimum internal shock radius, r ~ F^ro, is also the 
radius at which an individual sub-shell may experience significant change in its 
width ro, due to Lorentz factor variation across the shell. 

2.3 The allowed range of Lorentz factors and baryon loading 

The acceleration, 7 oc r, of fireball plasma is driven by radiation pressure. Fire- 
ball protons are accelerated through their coupling to the electrons, which are 
coupled to fireball photons. We have assumed in the analysis presented above, 
that photons and electrons are coupled throughout the acceleration phase. How- 
ever, if the baryon loading is too low, radiation may decouple from fireball elec- 
trons already at r < rj. The fireball Thomson optical depth is given by the 
product of comoving expansion time, r/j{r)c, and the photon Thomson scatter- 
ing rate, neCar- The electron and proton comoving number densities are equal, 
rie = rip, and are determined by equating the r independent mass flux carried 
by the wind, 4Trr^cj{r)npmp, to the mass loss rate from the underlying source, 
which is related to the rate L at which energy is emitted through M = L/^rjc^). 
Thus, during the acceleration phase, 7(r) = r/ro, the Thomson optical depth 
tt oc r~^. tt drops below unity at a radius r <rf = rjro if ry > ry*, where 

•'-(jjS^)" = >-»>'>°'^-V'. (") 

Here ro = 10''ro.7 cm. 

If 77 > 77* radiation decouples from the fireball plasma at 7 = r/ro = 
r]t^^ri~^^^. li rj ^ ry*, then most of the radiation energy is not converted to 
kinetic energy prior to radiation decoupling, and most of the fireball energy 
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escapes in the form of thermal radiation. Thus, the baryon load of fireball 
shells, and the corresponding final Lorentz factors, must be within the range 
10^ < r K -q < rjf: Ki \{)^ in order to allow the production of the observed 
non-thermal 7-ray spectrum. 



2.4 Fireball interaction with surrounding medium 

As the fireball expands, it drives a relativistic shock (blast-wave) into the sur- 
rounding gas, e.g. into the inter-stellar medium (ISM) gas if the explosion occurs 
within a galaxy. In what follows, we refer to the surrounding gas as "ISM gas," 
although the gas need not necessarily be inter-stellar. At early time, the fireball 
is little affected by the interaction with the ISM. At late time, most of the fire- 
ball energy is transferred to the ISM, and the flow approaches the self-similar 
blast- wave solution of Blandford & McKee |l5). At this stage a single shock 
propagates into the ISM, behind which the gas expands with Lorentz factor 

rsMir) = f TTT^) r-^" = 150 (^) r'^ . (12) 

where E = 10^^£'53 erg is the fireball energy, n = Iuq cm^^ is the ISM number 
density, and r = lO^'^rn cm is the shell radius. The characteristic time at which 
radiation emitted by shocked plasma at radius r is observed by a distant observer 



is t « r/4r|,,c |07| 



The transition to self-similar expansion occurs on a time scale T (measured 
in the observer frame) comparable to the longer of the two time scales set by the 
initial conditions: the (observer) GRB duration Tgrb and the (observer) time 
Tr at which the self-similar Lorentz factor equals the original ejecta Lorentz 
factor r, FsMit = Tp) = T. Since t = r/4r|4^c, 



T — max 



Tgrb, 5 ( r^ rj s 



(13) 



During the transition, plasma shocked by the reverse shocks expands with Lorentz 
factor close to that given by the self-similar solution, 

n,.^rBM{t = T)^2A5(-^\ T^^'\ (14) 

where T = lOTi s. The unshocked fireball ejecta propagate at the original expan- 
sion Lorentz factor, F , and the Lorentz factor of plasma shocked by the reverse 
shock in the rest frame of the unshocked ejecta is ~ r/FtY.. If T ~ Tgrb ^ Tr 
then F/FtY. ^ 1, the reverse shock is relativistic, and the Lorentz factor associ- 
ated with the random motion of protons in the reverse shock is 7^ ~ F/F^y . 

If, on the other hand, F Fr ^ Tgrb then F/Ft^. ~ 1, and the reverse 
shock is not relativistic. Nevertheless, the following argument suggests that the 
reverse shock speed is not far below c, and that the protons are therefore heated 
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to relativistic energy, 7^ — 1 ~ 1. The comoving time, measured in the fireball 
ejecta frame prior to deceleration, is ico. — r/Fc. The expansion Lorentz factor is 
expected to vary across the ejecta, AF/ F ^ 1, due to variability of the underlying 
GRB source over the duration of its energy release. Such variation would lead 
to expansion of the ejecta, in the comoving frame, at relativistic speed. Thus, 
at the deceleration radius, ico. — rT, the ejecta width exceeds ~ cico. — FcT. 
Since the reverse shock should cross the ejecta over a deceleration time scale, 
~ FT, the reverse shock speed must be close to c. We therefore conclude that the 
Lorentz factor associated with the random motion of protons in the reverse shock 
is approximately given by 7^ — 1 ~ F/Ft^. for both -T/Itr. ^ 1 and -T/Itr. S> 1. 

Since Tqrb 10 s is typically comparable to Tr, the reverse shocks are 
typically expected to be mildly relativistic. 



2.5 Fireball geometry 

We have assumed in the discussion so far that the fireball is spherically sym- 
metric. However, a jet-like fireball behaves as if it were a conical section of a 
spherical fireball as long as the jet opening angle is larger than F~^. This is 
due to the fact that the linear size of causally connected regions, cico. ^ r/F in 
the fireball frame, corresponds to an angular size ctco./r ^ F~^. Moreover, due 
to the relativistic beaming of radiation, a distant observer can not distinguish 
between a spherical fireball and a jet-like fireball, as long as the jet opening 
angle 9 > F^^. Thus, as long as we are discussing processes that occur when 
the wind is ultra-relativistic, 7^ ^ 300 (prior to significant fireball deceleration 
by the surrounding medium), our results apply for both a spherical and a jet- 
like fireball. In the latter case, L (E) in our equations should be understood as 
the luminosity (energy) the fireball would have carried had it been spherically 
symmetric. 



2.6 7-ray emission 

If the Lorentz factor variability within the wind is significant, internal shocks 
would reconvert a substantial part of the kinetic energy to internal energy. The 
internal energy may then be radiated as 7-rays by synchrotron and inverse- 
Compton emission of shock-accelerated electrons. The internal shocks are ex- 
pected to be "mildly" relativistic in the fireball rest frame, i.e. characterized by 
Lorentz factor I^i — 1 ~ a few. This is due to the fact that the allowed range 



of shell Lorentz factors is ~ 10 to ~ 10 (see §2^), implying that the Lorentz 
factors associated with the relative velocities are not very large. Since internal 
shocks are mildly relativistic, we expect results related to particle acceleration in 
sub-relativistic shocks (see [|4| for review) to be valid for acceleration in internal 
shocks. In particular, electrons are expected to be accelerated to a power law 
energy distribution, cine/d'ye oc for 7e > 7^, with p ~ 2 p 10 T^. 

The minimum Lorentz factor 7™ is determined by the following considera- 
tion. Protons are heated in internal shocks to random velocities (in the wind 
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frame) 7^ — — If electrons carry a fraction of the shock inter- 

nal energy, then 7^ ~ ^e(rnp/me). The characteristic frequency of synchrotron 
emission is determined by and by the strength of the magnetic field. As- 
suming that a fraction of the internal energy is carried by the magnetic field, 
4TTrfcr^B^/8TT = ^sLint., the characteristic observed energy of synchrotron pho- 
tons, = rh'^^^eB/rrieC, is 

e^.^HTej'ji^MeY. (15) 

In deriving Eq. (|l^) we have assumed that the wind luminosity carried by internal 
plasma energy. Lint., is related to the observed 7-ray luminosity through Lint. = 
L^/^e. This assumption is justified since the electron synchrotron cooling time 
is short compared to the wind expansion time (unless the equipartition fraction 
is many orders of magnitude smaller than unity), and hence electrons lose 
all their energy radiatively. Fast electron cooling also results in a synchrotron 
spectrum dn^/de^ oc e-y^ ^^"^ — at > e^i,, consistent with observed GRB 
spectra 

At present, there is no theory that allows the determination of the values 
of the equipartition fractions and £,b- Eq- (^H) implies that fractions not 
far below unity are required to account for the observed 7-ray emission. We 
note, that build up of magnetic field to near equipartition by electro-magnetic 
instabilities is expected to be a generic characteristic of collisionless shocks (see 
discussion in ref. and references therein), and is inferred to occur in other 
systems, e.g. in supernova remnant shocks (e.g. ||50| , p0| ). 

The 7-ray break energy e-^b of most GRBs observed by BATSE detectors is 



in the range of 100 keV to 300 keV [Ig]. It may appear from Eq. (15) that the 
clustering of break energies in this narrow energy range requires fine tuning of 
fireball model parameters, which should naturally produce a much wider range 
of break energies. This is, however, not the case |Q. Consider the dependence of 
e~fb on r. The strong F dependence of the pair-production optical depth, Eq. (^), 
implies that if the value of F is smaller than the minimum value allowed by Eq. 
(1), for which r^j^e-y = lOOMeV) « 1, most of the high energy photons in the 
power-law distribution produced by synchrotron emission, drij/de-y cx e^'^, would 
be converted to pairs. This would lead to high optical depth due to Thomson 
scattering on e^, and hence to strong suppression of the emitted flux p^ . For 
fireball parameters such that r-^^yie-y — lOOMeV) « 1, the break energy implied 



by Eqs. I\l§ and (|5|) is 

1/6 



Cyb^iefej'-^.MeY. (16) 



As explained in §2.3, shell Lorentz factors can not exceed 77* ~ 10'^, for which 
break energies in the X-ray range, e^b ~ 10 keV, may be obtained. We note, 
however, that the radiative flux would be strongly suppressed in this case too 
[p4|. If the typical F of radiation emitting shells is close to 77*, then the range 
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of Lorcntz factors of wind shells is narrow, which implies that only a small 
fraction of wind kinetic energy would be converted to internal energy which can 
be radiated from the fireball. 

Thus, the clustering of break energies at 1 MeV is naturally accounted for, 
provided that the variability time scale satisfies At < 10~^ s, which implies an 
upper limit on the source size, since At > tq/c. We note, that a large fraction of 
bursts detected by BATSE show variability on the shortest resolved time scale. 



10 ms 1 117]. In addition, a natural consequence of the model is the existence 



of low luminosity bursts with low, 1 keV to 10 keV, break energies |4j. Such 
"X-ray bursts" may have recently been identified . 

For internal collisions, the observed 7-ray variability time, ri/F'^c w At, 
reflects the variability time of the underlying source, and the GRB duration 
reflects the duration over which energy is emitted from the source. Since the 
wind Lorentz factor is expected to fluctuate on time scales ranging from the 
shortest variability time tq/c to the wind duration Tgrb , internal collisions will 
take place over a range of radii, r F'^tq to r ~ P^cTgrb- 



2.7 Afterglow emission 

Let us consider the radiation emitted from the reverse shocks, during the tran- 
sition to self-similar expansion. The characteristic electron Lorentz factor (in 
the plasma rest frame) is 7^ — (,e{r j rtr.)mp/ nie, where the internal energy 
per proton in the shocked ejecta is ~ {F / Ftr .)mpc'^ . The energy density U is 
E ss Anr'^cFF^j. U, and the number of radiating electrons is Ng « E/FnipC^. Us- 
ing Eq. (nj) and r — 4F^^ cF, the characteristic (or peak) energy of synchrotron 



photons (in the observer frame) is |113| 



e^^ « mr.ll^^ = 2ee,-i^]I'-inJ'rl, eV, (17) 
and the specific luminosity, — dL/de^, at is 

L„ « (2.a)-ir,,.-07Ve « lO'%%EU\-'/'F,:y,/'s-\ (18) 

where = 0.1^e,-i, and = 0.1^b,-i- 

Here too, we expect a power law energy distribution, dNe/d'ye cx 7"^ for 
7e > 7m, with p ~ 2. Since the radiative cooling time of electrons in the reverse 
shock is long compared to the ejecta expansion time, the specific luminosity 
extends in this case to energy > e^m as = Lm{e-y/ejm)~^^'^ , up to photon 
energy e^c- Here e-yc is the characteristic synchrotron frequency of electrons for 
which the synchrotron cooling time, dTririec/aTjeB^ , is comparable to the ejecta 
(rest frame) expansion time, ~ /tr.T'- At energy e-y > e^c, 

e,c « 0.^b'-We^3^'t-'/' keV, (19) 



the spectrum steepens to oc . 
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The shock driven into the ISM continuously heats new gas, and produces 
relativistic electrons that may produce the delayed afterglow radiation observed 
on time scales t 3> T, typically of order days to months. As the shock- wave 
decelerates, the emission shifts to lower frequency with time. Since we are inter- 
ested in proton acceleration to high energy and in the production of high energy 
neutrinos, which take place primarily in the internal and reverse shocks (see §^ 
) , we do not discuss in detail the theory of late-time afterglow emission. 



3 Some implications of afterglow observations 



Afterglow observations lead to the confirmation, as mentione d in the Introduc- 
tion, of the cosmological origin of GRBs j6l), and confirmed [106, 116 1 standard 
model predictions |^,^4|j72 , IOC | of afterglow that results from synchrotron emis- 
sion of electrons accelerated to high energy in the highly relativistic shock driven 
by the fireball into its surrounding gas. Since we are interested mainly in the 
earlier, internal collision phase of fireball evolution, we do not discuss afterglow 
observations in detail. We note, however, several implications of afterglow ob- 
servations which are of importance for the discussion of UHECR production. 

The following point should be clarified in the context of afterglow obser- 
vations. The distribution of GRB durations is bimodal, with broad peaks at 
Tgrb — 0.2 s and Tgrb = 20 s Q. The majority of bursts belong to the long 
duration, Tqrb ~ 20 s, class. The detection of afterglow emission was made pos- 
sible thanks to the accurate GRB positions provided on hour time scale by the 
BeppoSAX satellite [Q. Since the detectors on board this satellite trigger only 
on long bursts, afterglow observations are not available for the sub-population 
of short, Tgrb ~ 0.2 s, bursts. Thus, while the discussion of the fireball model 
presented in based on 7-ray observations and on simple phenomenological 
arguments, applies to both long and short duration bursts, the discussion below 
of afterglow observations applies to long duration bursts only. It should therefore 
be kept in mind that short duration bursts may constitute a different class of 
GRBs, which, for example, may be produced by a different class of progenitors 
and may have a different redshift distribution than the long duration bursts. 

Prior to the detection of afterglows, it was commonly assumed that the far- 
thest observed GRBs lie at redshift z ^ 1. Following the detection of afterglows 
and the determination of GRB redshifts, it is now clear that most GRB sources 
lie within the redshift range 2^0.5toz'~2, with some bursts observed at 
z > 3. For the average GRB 7-ray fluence, 1.2 x 10~^erg/cm^ in the 20 keV to 
2 MeV band, this implies characteristic isotropic 7-ray energy and luminosity 

~ 10^3 erg and ~ 10^^ gj.g/s (in the 20 keV to 2 MeV band), about an 
order of magnitude higher than the values assumed prior to afterglow detec- 
tion (Here, and throughout the paper we assume a flat universe with f2 — 0.3, 
A = 0.7, and Hq = 65km/sMpc). These estimates are consistent with more de- 
tailed analyses of the GRB luminosity function and redshift distribution. Mao 
& Mo, e.g., find, for the cosmological parameters we use, a median GRB energy 
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of « 0.6 X 10'^'^ erg in the 50keV to 300keV band |Q, corresponding to a median 
GRB energy of « 2 x lO^^erg in the 20 keV to 2 MeV band. 

The determination of GRB redshifts also lead to a modification of GRB rate 
estimates. Since most observed GRB sources lie within the redshift range z ^ 0.5 
to z 2, observations essentially determine the GRB rate per unit volume at 
z ~ 1. The observed rate of 10^ /yr imphes Rgrb{z = 1) « S/Gpc'^yr. The 
present, z = 0, rate is less well constrained, since available data are consistent 
with both no evolution of GRB rate with redshift, and with strong evolution 
(following, e.g., star formation rate), in which Rgrb{z = 1)/Rgrb{z = 0) ~ 10 
[ pOUSSl] . Detailed analyses, assuming i?GRB is proportional to star formation rate, 
lead to Rgrb{z = 0) ~ 0.5/Gpc'^yr ||6^. The implied local (z — 0) 7-ray energy 
generation rate by GRBs in the 20 keV to 2 MeV band is therefore 

e^(z = 0) = lO'^'^Cerg/MpcV, (20) 

with C in the range of ^ 10^°'^ to ^ 10°'^. Note, that is independent of the 
fireball geometry. If fireballs are conical jets of solid angle AH, then the total 
energy released by each burst is smaller by a factor Af^/An than the isotropic 
energy, and the GRB rate is larger by the same factor. 

Due to present technical limitations of the experiments, afterglow radiation 
is observed in most cases only on time scale >> 10 s. At this stage, radiation is 
produced by the external shock driven into the surrounding gas, and afterglow 
observations therefore do not provide direct constraints on plasma parameters 
at the internal and reverse shocks, where protons are accelerated to ultra-high 
energy. In one case, however, that of GRB 990123, optical emission has been 
detected on ~ 10 s time scale [Q. The most natural explanation of the observed 
optical radiation is synchrotron emission from electrons accelerated to high en- 
ergy in the reverse shocks driven into fireball ejecta at the onset of interaction 



with the surrounding medium |73,8^, as explained in §2.7. This observation pro- 
vides therefore direct constraints on the fireball ejecta plasma. First, it provides 
strong support for one of the underlying assumptions of the dissipative fireball 



scenario described in 1 2J1 , that the energy is carried from the underlying source 
in the form of proton kinetic energy. This is due to the fact that the observed 
radiation is well accounted for in a model where a shock propagates into fire- 
ball plasma composed of protons and electrons (rather than, e.g., pair plasma). 
Second, comparison of the observed flux with model predictions, Eqs. (|l^) and 
(^, implies ~ 10-1. 

Afterglow observations imply that a significant fraction of the energy initially 
carried by the fireball is converted into 7-rays, i.e. that the observed 7-ray energy 
provides a rough estimate of the total fireball energy. This has been demonstrated 
for one case, that of GRB970508, by a comparison of the total fireball energy 
derived from long term radio observations with the energy emitted in 7-rays 



[114 3J], and for a large number of bursts by a comparison of observed 7-ray 
energy with the total fireball energy estimate based on X-ray afterglow data |Q . 
In the context of the fireball model described in the inferred high radiative 
efficiency implies that a significant fraction of the wind kinetic energy must be 
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converted to internal energy in internal shocks, and that electrons must carry a 
significant fraction of the internal energy, i.e. that should be close to unity. 
We have already shown, see §2.6, that values not far below unity are required 
to account for the observed 7-ray emission. Conversion in internal shocks of a 
large fraction of fireball kinetic energy to internal energy is possibly provided 
the variance in the Lorentz factors of fireball shells is large |Q . 

In accordance with the implications of afterglow observations, we assume in 
the discussion below of UHECR and neutrino production in GRB fireballs, that 
the fraction of fireball energy converted to internal energy carried by electrons, 
and hence to 7-rays, is large. For discussion of UHECR production under the 
assumption that only a negligible fraction, ~ m^/mp, of the fireball energy is 
converted to radiation see ref. BSl. 



4 UHECRs from GRB fireballs 



4.1 Fermi acceleration in GRBs 



In the fireball model, the observed radiation is produced, both during the GRB 
and the afterglow, by synchrotron emission of shock accelerated electrons. In 
the region where electrons are accelerated, protons are also expected to be shock 
accelerated. This is similar to what is thought to occur in supernovae remnant 
shocks, where synchrotron radiation of accelerated electrons is the likely source of 
non-thermal X-rays (recent ASCA observations give evidence for acceleration of 
electrons in the remnant of SN1006 to lO-'^'^eV [Q), and where shock acceleration 
of protons is believed to produce cosmic rays with energy extending to ^ lO^^eV 
(see, e.g., for review). Thus, it is likely that protons, as well as electrons, are 
accelerated to high energy within GRB fireballs. Let us consider the constraints 
that should be satisfied by the fireball parameters in order to allow acceleration 
of protons to ^ 10^" eV. 

We consider proton Fermi acceleration in fireball internal shocks, which take 
place as the fireball expands over a range of radii, r ~ F^r^ to r ^ /^^cTgrb, 
and at the reverse shocks driven into fireball ejecta due to interaction with sur- 



rounding medium at r F^cT ^ F^cFqyib (see §2^, §2^). Both internal and 
reverse shocks are, in the wind rest-frame, mildly relativistic, i.e. characterized 
by Lorentz factors Fi — l ^ 1. Moreover, since reverse shocks do not cause strong 



deceleration of fireball plasma, see §2.4, the expansion Lorentz factor Ptr. of fire- 
ball plasma shocked by reverse shocks is similar to the fireball Lorentz factor F 
prior to interaction with the surrounding medium. Thus, plasma parameters, e.g. 
energy and number density, in the reverse shocks are similar to those obtained 
in internal shocks due to variability on time scale At ~ T. Results obtained 
below for internal shocks are therefore valid also for reverse shocks, provided At 
is replaced with T. 

Since the shocks we are interested in are mildly relativistic, we expect results 
related to particle acceleration in sub-relativistic shocks (see [Q for review) 
to be valid for our scenario. The predicted energy distribution of accelerated 
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protons is therefore dup/dtp oc e^^ similar to the predicted electron 

energy spectrum, which is consistent with the observed photon spectrum (see 

The most restrictive requirement, which rules out the possibility of acceler- 
ating particles to energy ~ 10^° eV in most astrophysical objects, is that the 
particle Larmor radius Rl should be smaller than the system size [Q. In our 
scenario we must apply a more stringent requirement, namely that should 
be smaller than the largest scale / over which the magnetic field fluctuates, since 
otherwise Fermi acceleration may not be efficient. We may estimate I as follows. 
The comoving time, i.e. the time measured in the wind rest frame, is t — r/Fc. 
Thus, regions separated by a comoving distance larger than r/F are causally 
disconnected, and the wind properties fluctuate over comoving length scales up 
to I ^ r/F. We must therefore require Rl < r/F. A somewhat more stringent re- 
quirement is related to the wind expansion. Due to expansion the internal energy 
is decreasing and therefore available for proton acceleration (as well as for 7-ray 
production) only over a comoving time t ^ r/Fc. The typical Fermi acceleration 
time is ta = /Rl/cP'^ | |5^ , [T^ , where /3c is the Alfven velocity and / '--^ 1 In 
our scenario /? ~ 1 leading to the requirement /Rl < r/F. This condition sets a 
lower limit to the required comoving magnetic field strength. Using the relations 
Rl = s'p/eB = ep/ FeB, where e'p = ep/ F is the proton energy measured in the 
fireball frame, and inr^ cF^ / 8tt — S^sLy/S^e, the constraint fR^ < r/F may 
be written as [ 

> 0.02rFl,el,,L-',„ (21) 

where Cp = 10^°ep^20 eV is the accelerated proton energy. Note, that this con- 
straint is independent of r, the internal collision radius. 

The accelerated proton energy is also limited by energy loss due to syn- 
chrotron radiation and interaction with fireball photons. As discussed in §6, the 
dominant energy loss process is synchrotron cooling. The condition that the 
synchrotron loss time, tgy = {dirmpC^ /aTm1)ep^B~^ , should be smaller than 
the acceleration time sets an upper limit to the magnetic field strength. Since 
the equipartition field decreases with radius, B^.p. oc r~^, the upper limit on 
the magnetic field may be satisfied simultaneously with (pl|) provided that the 
internal collisions occur at large enough radius 

r > 10''fF,~iel ,,cm. (22) 

Since collisions occur at radius r « F^cAt, the condition ( p^ is equivalent to a 
lower limit on F 

F > liQf''\l%At-T. (23) 

From Eqs. ( ^l| ) and (p3|), we infer that a dissipative ultra- relativistic wind, 
with luminosity and variability time implied by GRB observations, satisfies the 
constraints necessary to allow the acceleration of protons to energy > 10^° eV, 
provided that the wind bulk Lorentz factor is large enough, F > 100, and that 
the magnetic field is close to equipartition with electrons. The former constraint, 
F > 100, is remarkably similar to that inferred based on the 7-ray spectrum. 



104 1, 
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and r ^ 300 is the "canonical" value assumed in the fireball model. The latter 
constraint, magnetic field close to equipartition, must be satisfied to account for 
both 7-ray emission (see § ^.6D and afterglow observations (see §||). 

Finally, two points should be clarified. First, it has recently been claimed 
that ultra-high energy protons would lose most of their energy adiabatically, i.e. 
due to expansion, before they escape the fireball |Q . This claim is based on the 
assumptions that internal shocks, and therefore proton acceleration, occur at 
r ~ r^ro only, and that subsequently the fireball expands adiabatically. Under 
these assumptions, protons would lose most their energy by the time they escape. 



However, as emphasized both in this section and in § p.2| , internal shocks are 
expected to occur over a wide range of radii, and in particular at r F^cT during 
the transition to self-similar expansion. Thus, proton acceleration to ultra-high 
energy is expected to operate over a wide range of radii, from r ^ i^^^o up to 
r ~ r^cT, where ultra-high energy particles escape. 

Second, it has recently been claimed in ji^ that the conditions at the external 
shock driven by the fireball into the ambient gas are not likely to allow proton 
acceleration to ultra-high energy. Regardless of the validity of this claim, it is 
irrelevant for the acceleration in internal shocks, the scenario considered for 
UHECR production in GRBs in both |104| and Q. Moreover, it is not at 
all clear that UHECRs can not be produced at the external shock, since the 
magnetic field may be amplified ahead of the shock by the streaming of high 
energy particles. For discussion of high energy proton production in the external 
shock and its possible implications see ref. |29] . 



4.2 UHECR flux and spectrum 



119 93] results confirm the flattening of the 



31 and AGASA \ 

10^^ eV, evidence for which existed in previous ex- 



Fly's Eye |1 
cosmic-ray spectrum at 
periments with weaker statistics |102]. Fly's Eye data is well fitted in the energy 
range 10^^-^ eV to 10^^'^ eV by a sum of two power laws: A steeper component, 
with differential number spectrum J oc E^^'^'^, dominating at lower energy, and 
a shallower component, J oc E~'^-^^, dominating at higher energy, E > 10^^ eV. 
The flattening of the spectrum, combined with the lack of anisotropy and the 
evidence for a change in composition from heavy nuclei at low energy to light 
nuclei (protons) at high energy 1 102 , 12|jl^,^,^ , suggest that a Galactic compo- 
nent of heavy nuclei, J cx E~'^^^'^, dominates the cosmic-ray flux at low energy, 
while an extra-Galactic component of protons, J oc E^'^-^^, dominates the flux 
at high energy, > 10^® eV. 

The GRB energy observed in 7-rays reflects the fireball energy in accelerated 
electrons. If accelerated electrons and protons carry similar energy, as indicated 
by afterglow observations |Q (see then the GRB cosmic-ray production 
rate is [see Eq. (|20|)] 



el{dnp/dep)2=o » lO'^'^erg/Mpc yr. 



(24) 
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Fig. 1. The UHECR flux expected in a cosmological model, where high-energy protons 
are produced at a rate {e^dhp / de^) = 0.8x lO^'^erg/Mpc^yr as predicted by the GRB 
model [Eq. @], solid line, compared to the Fly's Eye jl3|, Yakutsk Q and AGASA 
]4e| ] data, la flux error bars are shown. The highest energy points are derived assuming 
the detected events represent a uniform flux over the energy range 10^° eV-3 X 10^° eV. 
The dashed line is the sum of the GRB model flux and the Fly's Eye fit to the Galactic 
heavy nuclei component, J oc (~^'^ [Il3| (with normalization increased by 25%). 



In Fig. 1 we compare the UHECR spectrum, reported by the Fly's Eye p^ , 
the Yakutsk |Q, and the AGASA experiments with that expected from a 
homogeneous cosmological distribution of sources, each generating a power law 
differential spectrum of high energy protons dn/dep oc e~^. The absolute flux 
measured at 3 x lO^'* cV differs between the various experiments, corresponding 
to a systematic ~ 10% (~ 20%) over-estimate of event energies in the A GAS A 
(Yakutsk) experiment compared to the Fly's Eye experiment (see also |ll9|] ). 
In Fig. 1, the Yakutsk energy normalization is used. For the model calculation, 
a flat universe, SI = 0.3, A — 0.7 and Hp — 65km/Mpc s were assumed. The 
calculation is similar to that described in |l05|] . The generation rate of cosmic- 
rays (per unit comoving volume) was assumed to evolve rapidly with redshift 
following the luminosity density evolution of QSOs ||l7||, which is also similar to 
that describing the evolution of star formation rate [|64 67 : ncR{z) oc (1 -|- z)" 
with a w 3 1^ at low redshift, z < 1.9, ncR.{z) — Const, for 1.9 < z < 2.7, 
and an exponential decay at z > 2.7 |90). The cosmic-ray spectrum at energy 
> lO""^^ eV is little affected by modifications of the cosmological parameters or of 
the redshift evolution of cosmic-ray generation rate. This is due to the fact that 
cosmic-rays at this energy originate from distances shorter than several hundred 
Mpc. The spectrum and flux at Cp > 10^^ eV is mainly determined by the present 
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(z = 0) generation rate and spectrum, which in the model shown in Fig. 1 is 
elidhp/de),^a = 0.8 x lO^^g^g/Mpc^r- 

The suppression of model flux above 10^^ '' eV is due to energy loss of high 
energy protons in interaction with the microwave background, i.e. to the "GZK 
cutoff" ^,12C|. The available data do not allow to determine the existence (or 



absence) of the "cutoff" with high confidence. The AGASA results show an 
excess (at a ^ 2.5a confidence level) of events compared to model predictions 
above 10^°eV. This excess is not confirmed, however, by the other experiments. 
Moreover, since the 10^'^eV flux is dominated by sources at distances < 100 Mpc, 
over which the distribution of known astrophysical systems (e.g. galaxies, clus- 
ters of galaxies) is inhomogeneous, significant deviations from model predictions 
prese nted in Fig. 1 for a uniform source distribution are expected at this energy 
[ |105| ■ Clustering of cosmic-ray sources leads to a standard deviation, tr, in the 
expected number, N, of events above 10^° eV, given by a/N = 0.9(do/10Mpc)°-^ 
[0 , where do is the unknown scale length of the source correlation function and 
do ^ 10 Mpc for field galaxies. 

Thus, GRB fireballs would produce UHECR flux and spectrum consistent 
with that observed, provided the efficiency with which the wind kinetic energy is 
converted to 7-rays, and therefore to electron energy, is similar to the efficiency 



with which it is converted to proton energy, i.e. to UHECRs |104|. There is 



however, one additional point which requires consideration |104 . The energy of 
the most energetic cosmic ray detected by the Fly's Eye experiment is in excess 
of 2 X lO^^eV, and that of the most energetic AGASA event is ~ 2 x 10^°eV. 
On a cosmological scale, the distance traveled by such energetic particles is 
smafi: < lOOMpc (50Mpc) for the AGASA (Fly's Eye) event (e.g., 0). Thus, 
the detection of these events over a ~ Syr period can be reconciled with the rate 
of nearby GRBs, ^ 1 per 100 yr to ~ 1 per 1000 yr out to lOOMpc, only if there 
is a large dispersion, > lOOyr, in the arrival time of protons produced in a single 
burst (This implies that if a direct correlation between high energy CR events 
and GRBs, as suggested in is observed on a ~ lOyr time scale, it would be 
strong evidence against a cosmological GRB origin of UHECRs). 

The required dispersion is likely to occur due to the combined effects of de- 



flection by random magnetic fields and energy dispersion of the particles | 104 |. 
Consider a proton of energy propagating through a magnetic field of strength 
B and correlation length A. As it travels a distance A, the proton is typically de- 
flected by an angle a ~ X/Rl, where Rl = ep/eB is the Larmor radius. The typ- 
ical deflection angle for propagation over a distance D is 9s ~ {2D/9Xy^^ X/ Rl- 
This deflection results in a time delay, compared to propagation along a straight 
line, 

T(ep, D) « elD/4c « 10^e;;i,Dl,,XMr.cBU yr, (25) 

where D = lOOI^iooMpc, A = IAmpc Mpc and B = 10"^B_8 G. Here, we 
have chosen numerical values corresponding to the current upper bound on the 
inter-galactic magnetic field, BA^/^ < 10"^G Mpc^^^ 111,111 • The upper bound 
on the (systematic increase with redshift of the) Faraday rotation measure of 
distant, z < 2.5, radio sources, RM < 5rad/m^, implies an upper bound B < 
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10^^^(/i/0.75)(i7f,/i^)^^ G on an inter-galactic field coherent over cosmological 
scales [Q. Here, h is the Hubble constant in units of lOOkm/sMpc and fib is 
the baryon density in units of the closure density. For a magnetic field coherent 
on scales - A, this implies SA^/^ < lO-S{h/Omy/^{f2bh^ /OM)-^G Mpc^^^ 

The random energy loss UHECRs suffer as they propagate, owing to the 
production of pions, implies that at any distance from the observer there is 
some finite spread in the energies of UHECRs that are observed with a given 
fixed energy. For protons with energies > lO'^^eV the fractional RMS energy 
spread is of order unity over propagation distances in the range 10 — lOOMpc 
(e.g. |l|). Since the time delay is sensitive to the particle energy, this implies that 
the spread in arrival time of UHECRs with given observed energy is comparable 
to the average time delay at t hat e nergy T(ep, D) (This result has been confirmed 
by numerical calculations in [112|). Thus, the required time spread, r > 100 yr, 
is consistent with the upper bound, r < 10^ yr, implied by the present upper 
bound to the inter-galactic magnetic field. 



5 GRB model predictions for UHECR experiments 
5.1 The Number and Spectra of Bright Sources 

The initial proton energy, necessary to have an observed energy Cp, increases with 
source distance due to propagation energy losses. The rapid increase of the initial 
energy after it exceeds, due to electron-positron production, the threshold for 
pion production effectively introduces a cutoff distance, Dc{ep), beyond which 
sources do not contribute to the fiux above ep. The function Dc{ep) is shown 
in Fig. H (adapted from Since Dc{ep) is a decreasing function of ep, for a 

given number density of sources there is a critical energy Ec, above which only 
one source (on average) contributes to the flux. In the GRB model tc depends on 
the product of the burst rate Rgrb and the time delay. The number of sources 



contributing, on average, to the flux at energy tp is [76 

Ait 



N{ep) = —RGRBD,{epYT[ep,D,{ep)] , (26) 



and the average intensity resulting from all sources is 



J{ep) ^ ^RGRB^^D,{ep) , (27) 

where dup/dep is the number per unit energy of protons produced on average by 
a single burst (this is the formal definition of Dc{ep)). The critical energy ec is 
given by 

A'jr 

—RGRBD,{e,)\[e,,D,{e,)]^l . (28) 
5 

Cc, the energy beyond which a single source contributes on average to the 
flux, depends on the unknown properties of the inter-galactic magnetic field, 
T (X B^\. However, the rapid decrease of Dc{ep) with energy near 10^°eV implies 
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Fig. 2. Ec, the energy beyond which a single GRB contributes on average to the UHECR 
flux, as a function of the product of GRB rate, -Rgrb ~ 1/Gpc^, and the time delay 
of a lO'^" eV proton originating at 100 Mpc distance. The time delay depends on the 
unknown inter-galactic field, t oc B^X. Dashed lines show the allowed range ofB^A: The 
lower limit is set by the requirement that at least a few GRB sources be present at D < 
100 Mpc, and the upper limit by the Faraday rotation bound BX^^^ < 10"**G Mpc^^^ 
|5q], see Eq. (Gq) and the discussion the follows it. 



see 
„ is 



that Ec is only weakly dependent on the value of B^X, as shown in Fig. ^. In The 
GRB model, the product Rgrbt{D — lOOMpc, Cp = lO^^eV) is approximately 
limited to the range 10~^ Mpc~'^ to 10~^ Mpc~^ [The lower limit is set by the 
requirement that at least a few GRB sources be present at D < 100 Mpc, and 
the upper limit by the Faraday rotation bound BX^^^ < 10~*G Mpc^^^ |5^, 
Eq. (|2^) , and Rgrb ^1/ Gpc'^yr] . The corresponding range of values of e, 
lO^OeV < ec < 4 X lO^OeV. 

Fig. ^ presents the flux obtained in one realization of a Monte-Carlo simula- 
tion described by Miralda-Escude & Waxman of the total number of UHE- 
CRs received from GRBs at some fixed time. For each realization the positions 
(distances from Earth) and times at which cosmological GRBs occurred were ran- 
domly drawn, assuming an intrinsic proton generation spectrum dup/dep oc e" ^, 
and Ec = 1-4 X lO^^eV. Most of the realizations gave an overall spectrum similar 
to that obtained in the realization of Fig. || when the brightest source of this real- 
ization (dominating at 10^'^eV) is not included. At < ec, the number of sources 
contributing to the flux is very large, and the overall UHECR flux received at 
any given time is near the average (the average flux is that obtained when the 
UHECR emissivity is spatially uniform and time independent). At ep > ec, the 
flux will generally be much lower than the average, because there will be no 
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Fig. 3. Results of a Monte-Carlo realization of the bursting sources model, with = 
1.4 X lO'^" eV: Thick solid line- overall spectrum in the realization; Thin solid line- 
average spectrum, this curve also gives Dc(ep); Dotted lines- spectra of brightest sources 
at different energies. 



burst within a distance Dc{ep) having taken place sufRciently recently. There is, 
however, a significant probabihty to observe one source with a flux higher than 
the average. A source similar to the brightest one in Fig. || appears ~ 5% of the 
time. 

At any fixed time a given burst is observed in UHECRs only over a narrow 
range of energy, because if a burst is currently observed at some energy ep then 
UHECRs of much lower energy from this burst have not yet arrived, while higher 
energy UHECRs reached us mostly in the past. As mentioned above, for energies 
above the pion production threshold, ep ~ 5 x lO^^eV, the dispersion in arrival 
times of UHECRs with fixed observed energy is comparable to the average delay 
at that energy. This implies that the spectral width Aep of the source at a given 
time is of order the average observed energy, Aep ~ ep. Thus, bursting UHECR 
sources should have narrowly peaked energy spectra, and the brightest sources 
should be different at different energies. For steady state sources, on the other 
hand, the brightest source at high energies should also be the brightest one at low 
energies, its fractional contribution to the overall flux decreasing to low energy 
only as Dc{ep)^^. A detailed numerical analysis of the time dependent energy 
spectrum of bursting sources is given in | ]9^ , p3| . 
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5.2 Spectra of Sources at Cp < 4 X lO^^eV 

The detection of UHECRs above lO^^eV imply that the brightest sources must 
he at distances smaller than lOOMpc. UHECRs with ep < 4 x lO^^eV from 
such bright sources will suffer energy loss only by pair production, because at 
ep < 5 X 10^^ eV the mean-free-path for pion production interaction (in which 
the fractional energy loss is ~ 10%) is larger than IGpc. Furthermore, the energy 
loss due to pair production over lOOMpc propagation is only ~ 5%. 

In the case where the typical displacement of the UHECRs due to deflections 
by inter-galactic magnetic fields is much smaller than the correlation length, 
A > DOsiD, Cp) ~ D{D/X)^/^X/Rl, all the UHECRs that arrive at the observer 
are essentially deflected by the same magnetic field structures, and the absence 
of random energy loss during propagation implies that all rays with a fixed 
observed energy would reach the observer with exactly the same direction and 
time delay. At a fixed time, therefore, the source would appear mono-energetic 
and point-like. In reality, energy loss due to pair production results in a finite 
but small spectral and angular width, Atp/ tp ' 



< 1% [115 



In the case where the typical displacement of the UHECRs is much larger 
than the correlation length, A ^ D9s{D, ep), the deflection of different UHECRs 
arriving at the observer are essentially independent. Even in the absence of any 
energy loss there are many paths from the source to the observer for UHECRs 
of fixed energy that are emitted from the source at an angle < Og relative to 
the source-observer line of sight. Along each of the paths, UHECRs are defiected 
by independent magnetic field structures. Thus, the source angular size would 
be of order 9s and the spread in arrival times would be comparable to the 
characteristic delay r, leading to Atp/ tp ^ 1 even when there are no random 
energy losses. The observed spectral shape of a nearby {D < lOOMpc) bursting 
source of UHECRs at < 4 x lO^^eV was derived for the case A <^ D9s{D, ep) 



in 1 115 1, and is given by 



dN 

' r» — 1 



cxp 



]it,D) 



(29) 



where ea{t,D) = De{2B'^X/3ct)^/'^. For this spectrum, the ratio of the RMS 
UHECR energy spread to the average energy is 30% 

Fig. ^ shows the line 9sD = A in the B — \ plane, for a source at a distance 
D = 30Mpc observed at energy Cp ~ lO^^eV. Since the 9sD = A line divides the 
allowed region in the plane at A ~ IMpc, measuring the spectral width of bright 
sources would allow to determine if the field correlation length is much larger, 
much smaller, or comparable to IMpc. 



6 High energy Neutrinos 

6.1 Internal shock (GRB) neutrinos 

Neutrinos at energies ~ IQ-*-^ eV. Protons accelerated in the fireball to 
high energy lose energy through photo-meson interaction with fireball photons. 
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Fig. 4. The line 6sD = A for a source at 30Mpc distance observed at energy ~ 10^ eV 
(dot-dash line), shown with the Faraday rotation upper limit B\^''^ < 10"*G Mpc^^^ 
(sohd line), and with the lower limit BA^^^ > 10"^°G Mpc^''^ required in the GRB 
model [see Eq. (||)]. 



The decay of charged pions produced in this interaction, tt"*" —> fi'^ + Vfj. ^ 
e'^ + Ue +T'fj^ + v^, results in the production of high energy neutrinos. The key 
relation is between the observed photon energy, e-y, and the accelerated proton's 
energy, Cp, at the threshold of the Z\-resonance. In the observer frame, 

e-yep = 0.2 GeV^ . (30) 

For r « 300 and — 1 MeV, we see that characteristic proton energies ~ 
10^^ eV are required to produce pions. Since neutrinos produced by pion decay 
typically carry 5% of the proton energy (see below), production of ~ lO^"* eV 



neutrinos is expected |109,110 



The fractional energy loss rate of a proton with energy = ep/F measured 
in the wind rest frame due to pion production is 

i-i(,') 



1 /"^ /"CXD 

= TT"^ / dea^{e)^{e)e / dx x~'^ —^{e^ = x) , (31) 

where 7p = e'p/rapC^ , cr-n-Ce) is the cross section for pion production for a photon 
with energy e in the proton rest frame, ^(e) is the average fraction of energy lost 
to the pion, eq = 0.15GeV is the threshold energy, and dn^/de^ is the photon 
density per unit photon energy in the wind rest frame. In deriving Eq. (|3l|) we 
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have assumed that the photon distribution in the wind rest frame is isotropic. 
The GRB photon spectrum is well fitted in the BATSE detector range, 20 keV 
to 2 MeV, by a combination of two power-laws, dn^/dej oc e"'', with /3 ~ 1 at 
e-y < e-yb, /3 ~ 2 at > e-yb and e^b IMeV Q. Thus, the second integral in 



may be approximated by 



where Uj is the photon energy density (in the range corresponding to the ob- 
served BATSE range) in the wind rest-frame, /3 = 1 for e < e'^^ and /3 = 2 for 
e > d^^. is the break energy measured in the wind frame, d^^^ — e^b/T. 
The main contribution to the first integral in (|3l| ) is from photon energies 
e ^ Epoak — O.SGeV, where the cross section peaks due to the A resonance. 
Approximating the integral by the contribution from the resonance we obtain 

^^^(4) - ^TT-cCTpcak^pcak^^ min(l, 27pe!:j,^/epcak) • (33) 

^^-yb £pcak 

Here, cTpcak — 5 x 10~^^cm^ and ^pcak — 0.2 are the values of a and ^ at e = Cpeakj 
and Ae ~ 0.2GeV is the peak width. 

The energy loss of protons due to pion production is small during the accel- 
eration process. Once accelerated, the time available for proton energy loss by 
pion production is comparable to the wind expansion time as measured in the 
wind rest frame, t^o ~ r/Fc. Thus, the fraction of energy lost by protons to pions 
is /tt — r/FctTT- The energy density in the BATSE range, U.y, is related to the 
luminosity by L y = Anr^F-^cUj. Using this relation in (|3^), and r = 2F-^cAt, 
U is given by ||l09|| 

/^(ep) « 0.1- n4 A 4- ^ W' /<: . r^+W-„nJL' i^'^) 



e7b,McV^2'.5^i-2 I Ep/^pb' Otherwise. 
The proton break energy is 

epb ~ 10i6r2%(e^6,Mcv)"' cV . (35) 

The value of /^r, Eq. (|3^), is strongly dependent on F. It has recently been 
pointed out that if the Lorentz factor F varies significantly between bursts, 
with burst to burst variations AF/F ~ 1, then the resulting neutrino flux will 
be dominated by a few neutrino bright bursts, and may significantly exceed the 
flux implied by Eq. (|3^), derived for typical burst parameters. This may strongly 
enhance the detectability of G RB neutrinos by planned neutrino telescopes Q . 
However, as explained in § ^2.6^ the Lorentz factors of fireballs producing observed 
GRBs can not differ significantly from the minimum value allowed by Eq. (^), 
F ~ 250, for which the fireball pair production optical depth, Eq. (jj), is w 1 
for e.y = 100 MeV: Lower Lorentz factors lead to optically thick fireballs, while 
higher Lorentz factors lead to low luminosity X-ray bursts (which may have 
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already been identified). Thus, for Lorcntz factors consistent with observed GRB 
spectra, for which T^-y{e-y — lOOMeV) « 1), we find 

e^bMeV^t]/^ lep/epfc, otherwise. 

A detailed analysis, using Monte-Carlo simulations of the internal shock model, 
confirms that for fireball parameter range consistent with observed GRB char- 
acteristics, /tt at Cp > £pb is limited to the range of ~ 10% to 30% 

Thus, for parameters typical of a GRB producing wind, a significant fraction 
of the energy of protons accelerated to energies larger than the break energy, 
~ lO^^eV, would be lost to pion production. Roughly half of the energy lost 
by protons goes into tt" 's and the other half to 7r+ 's. Neutrinos are produced 
by the decay of tt+'s, 7r+ ^ ^+ + j/^ e+ + i^e + + t'/x [the large optical 
depth for high energy 7's from tt'~' decay, Eq. (||), would not allow these photons 
to escape the wind]. The mean pion energy is 20% of the energy of the proton 
producing the pion. This energy is roughly evenly distributed between the 7r+ 
decay products. Thus, approximately half the energy lost by protons of energy 
Cp is converted to neutrinos with energy ^ 0.05ep. Eq. (jsj) then implies that the 
spectrum of neutrinos above ej/f, = O.OSepb follows the proton spectrum, and is 
harder (by one power of the energy) at lower energy. 

If GRBs are the sources of UHECRS, then using Eq. (13) the expected GRB 



neutrino flux is |11C| 



c f 

« 1.5 X 10-9M|^min{l,e,/e,fc}GeVcm-2s-isr-i, (37) 

where tn ~ 10^° yr is the Hubble time. The factor of 1/8 is due to the fact that 
charged pions and neutral pions are produced with roughly equal probabilities 
(and each neutrino carries ~ 1/4 of the pion energy). 

The GRB neutrino flux can also be estimated directly from the observed 
gamma-ray fluence. The BATSE detectors measure the GRB fiuence over two 
decades of photon energy, ~ 0.02MeV to ~ 2MeV, corresponding to a decade 
of radiating electron energy (the electron synchrotron frequency is proportional 
to the square of the electron Lorentz factor). If electrons carry a fraction 
of the energy carried by protons, then the muon neutrino fluence of a single 
burst is eldN^/de„ « 0.25(/Tr/i^e)^7/ In(lO). The average neutrino flux per unit 
time and solid angle is obtained by multiplying the single burst fluence with 
the GRB rate per solid angle, w lO'^ bursts per year over 4tt sr. Using the 
average burst fluence — 10~^erg/cm^, we obtain a muon neutrino flux e^ip^ « 
3 X 10~^(/7r/Ce)GeV/cm^s sr. Thus, the neutrino flux estimated directly from 
the gamma-ray fluence agrees with the estimate ( ^7| ) based on the cosmic-ray 
production rate. 
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Neutrinos at energy > 10 eV. The neutrino spectrum (37) is modified 
at high energy, where neutrinos are produced by the decay of muons and pions 
whose hfe time t^,^ exceeds the characteristic time for energy foss due to adi- 



abatic expansion and synchrotron emission [109 8(:,110|. The synchrotron loss 
time is determined by the energy density of the magnetic field in the wind rest 
frame. For the characteristic parameters of a GRB wind, the muon energy for 
which the adiabatic energy loss time equals the muon life time, e^^, is comparable 
to the energy at which the life time equals the synchrotron loss time, r*. For 
pions, e° > e^. This, and the fact that the adiabatic loss time is independent 
of energy and the synchrotron loss time is inversely proportional to energy, im- 
ply that synchrotron losses are the dominant effect suppressing the flux at high 
energy. The energy above which synchrotron losses suppress the neutrino flux is 



^ 1/2 

^^(P^.l.,) ^ (£.B J \ rn2 A , , / 10, for V, 



e,, -^VCe^^-'V A.5^i-2e^^McV X I ^QQ^ ' (38) 

The efficiency of neutrino production in internal collisions decreases with At, 
U oc At-^ [see Eq. (||)], since the radiation energy density is lower at larger 
collision radii. However, at larger radii synchrotron losses cut off the spectrum 



at higher energy, e'^{At) oc At [see Eq. (38)]. Collisions at large radii therefore 
result in extension of the neutrino spectrum of Eq. ( |37| ) to higher energy, beyond 
the cutoff energy Eq. (^Sj), 

el<P^ oc e^>el. (39) 



Comparison with other authors. We note, that the results presented above 
were derived using the Z\-approximation, i.e. assuming that photo-meson interac- 
tions are dominated by the contribution of the Z\-resonance. It has recently been 
shown 1?^, that for photon spectra harder than drij/dej oc e~^, the contribution 
of non-resonant interactions may be important. Since in order to interact with 
the hard part of the photon spectrum, e-y < e^b, the proton energy must exceed 
the energy at which neutrinos of energy e^b are produced, significant modifica- 
tion of the Z\- approximation results is expected only for e,y ^ e^b, where the 
neutrino flux is strongly suppressed by synchrotron losses. 

The neutrino flux from GRBs is small above lO^^eV, and a neutrino flux 
comparable to the 7-ray flux is expected only below ~ lO^^eV, in agreement 
with the results of rcf. [Q. Our result is not in agreement, however, with that 
of ref. [101[, where a much higher flux at ~ lO^^eV is obtained based on the 
equations of ref. [109[, which are the same equations as used hereR. 



^ The parameters chos en i n |101| are Lj = 10^'^erg/s, At — 10s, and F — 100. Using 
equation (4) of ref. [109|, which is the same as Eq. (^) of the present paper, we 
obtain for these parameters = 1.6 x 10~*, while the author of [101| obtains, using 
the same equation, /,r = 0.03. 
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6.2 Reverse shock (afterglow) neutrinos, ~ lO-*^* eV 

Let us now consider neutrino emission from photo-meson interactions of pro- 
tons accelerated to high energies in the reverse shocks driven into the fireball 
ejecta at the initial stage of interaction of the fireball with its surrounding gas, 
which occurs on time scale T ^ 10 s, comparable to the duration of the GRB 
itself (see §2.4). Optical-UV photons are radiated by electrons accelerated in 
shocks propagating backward into the ejecta (see §2.7), and may interact with 
accelerated protons. The interaction of these low energy, 10 eV-1 keV, photons 
and high energy protons produces a burst of duration ~ T of ultra-high energy, 
10^^-10^^ eV, neutrinos [as indicated by Eq. (|30|)] via photo-meson interactions 

[In). 

Afterglows have been detected in several cases; reverse shock emission has 
only been identified for GRB 990123 Both the detections and the non- 
detections are consistent with shocks occurring with typical model parameters 
[ 113 89| , |73| , suggesting that reverse shock emission may be common. The pre- 
dicted neutrino emission depends, however, upon parameters of the surrounding 
medium that can only be estimated once more observations of the prompt opti- 
cal afterglow emission are available. We first consider the case where the density 
of gas surrounding the fireball is n ~ lcm~'^, a value typical to the inter-stellar 
medium and consistent with GRB 990123 observations. 

The photon density in Eq. ( |3l| ) is related to the observed specific luminosity 
by dn-y/de^{x) — L^l^Fx) j (Airr'^crx). For proton Lorentz factor €a/2e'^^ ^ 7p < 
eo/2e^^^, where primed energies denote rest frame energies (e.g. — e^m/ Ttr.), 
photo-meson production is dominated by interaction with photons in the energy 

— 1/2 I I 

range e^™ < ^ e-^c, where oc e-y (see §2.7). For this photon spectrum, 
the contribution to the first integral of Eq. (|l]) from photons at the Z\-resonance 
is comparable to that of photons of higher energy, and we obtain 

,—1/ /\ 2 / Lm f £pcak \ CTpcak^pcak'^e ( Ar\\ 

Ar. is the expansion Lorentz factor of plasma shocked by the reverse shocks, 
given by Eq. (|lj). The time available for proton energy loss by pion production is 
comparable to the expansion time as measured in the wind rest frame, '~ r/Ftr.c. 
Thus, the fraction of energy lost by protons to pions is 

« 0.1 (y^^) ' X (e,™,evep,2o)^/^ (41) 

Eq. (^Tj) is valid for protons in the energy range 

4 X 10i« (^^^ (6^,,kcv)-'eV < ep < 4 X lO^i (^^^ (6^„,^ev)~'eV . (42) 
Such protons interact with photons in the energy range e^m to e-yc, where the 

1/2 

photon spectrum oc e-yc and the number of photons above interaction thresh- 

1/2 

old is cx Cp . At lower energy, protons interact with photons of energy > e^c-, 
rather then oc e~^/^. At these energies therefore /^r oc ep. 
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Since approximately half the energy lost by protons of energy is con- 
verted to neutrinos with energy ^ O.OSep, Eq. ( ^2[ ) implies that the spectrum 
of neutrinos below e^h « 10^^(rtr./250)^(e^c,kov)^^cV is harder by one power 
of the energy then the proton spectrum, and by half a power of the energy 
at higher energy. For a power law differential spectrum of accelerated protons 
dup/dcp oc e~^, the differential neutrino spectrum is dn^/dci, oc e~" with a = 1 
below the break and a = 3/2 above the break. Assuming that GRBs are indeed 
the sources of ultra-high energy cosmic rays, then Eqs. (|4^,^) and (|2j) imply 
that the expected neutrino intensity is 

" ^^""^ (lO^fcv) GeVcm-s-sr-, (43) 

where /i^^' = A(ep,2o = 2) and /3 = 1/2 for > 10"eV and /? = 1 for 
< IQi^eV. 

Some GRBs may result from the collapse of a massive star, in which case the 



fireball is expected to expand into a pre-existing wind (e.g. [^,118|). For typical 
wind parameters, the transition to self-similar behavior takes place at a radius 
where the wind density is n « 10'*cm~'^ ^ lcm~^. The higher density implies 
a lower Lorenz factor of the expanding plasma during the transition stage, and 
a larger fraction of proton energy lost to pion production. Protons of energy 
Cp > 10^^ eV lose all their energy to pion production in this case. If most GRBs 
result from the collapse of massive stars, then the predicted neutrino flux is 
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ei-^^ « 10"**min{l,e°'' /10'W}GeVcm~''s"'sr"\ (44) 

The neutrino flux is expected to be strongly suppressed at energy ei, > 
10^^ eV, since protons are not expected to be accelerated to energy :s> 10^° eV. 
If protons are accelerated to much higher energy, the {vp^^ v^) flux may extend 
to ~ W^n^^'^Q^Ll eV (~ lO^On-'/'^-'/' eV). At higher energy, synchrotron 
losses of pions (muons) will suppress the neutrino flux. 



6.3 Inelastic p-n collisions 

The acc elera tion, 7 oc r, of fireball plasma emitted from the source of radius 



ro (see §2^) is driven by radiation pressure. Fireball protons are accelerated 
through their coupling to the electrons, which are coupled to fireball photons. 
Fireball neutrons, which are expected to exist in most progenitor scenarios, are 
coupled to protons by nuclear scattering as long as the comoving p-n scattering 
time is shorter than the comoving wind expansion time r/^c — r^/c. As the 
fireball plasma expands and accelerates, the proton density decreases, Up (X 
r~^7~^, and neutrons may become decoupled. For rj > r/pn, where 

Vpn « mLli\-y\ (45) 
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neutrons decouple from the accelerating plasma prior to saturation, "f — r], 
at r — 'qpl^T]^^^'^ In this case, relativistic relative velocities between 

protons and neutrons arise, which lead to pion production through inelastic 
nuclear collisions. Since decoupling occurs at a radius where the collision time 
is similar to wind expansion time, each n leads on average to one pair of vv. 
The typical comoving neutrino energy, ~ 50 MeV, implies an observed energy 
~ 10 GeV. A typical burst, E = 10^^ erg at z = 1, with significant neutron 
to proton ratio and r/ > 400 will therefore produce a fluence F{i'e + Vg) ~ 
0.5F(iy^ + v^,) - lO-^cm-^ of ~ 10 GeV neutrinos. 

Relativistic relative p-n velocities, leading to neutrino production through 
inelastic collisions, may also result from diffusion of neutrons between regions of 
the fireball wind with large difference in _r |74|. If, for example, plasma expanding 
with very high Lorentz factor, F > 100, is confined to a narrow jet surrounded 
by a slower, r' ^ 10 wind, internal collisions within the slower wind can heat 
neutrons to relativistic temperature, leading to significant diffusion of neutrons 
from the slower wind into the faster jet. Such process may operate for winds 
with rj < 400 as well as for ry > 400, and may lead, for certain (reasonable) 
wind parameter values, to ^ 10 GeV neutrino flux similar to that due to p-n 
decoupling in a 77 > 400 wind. 



6.4 Implications 

The high energy neutrinos predicted in the dissipative wind model of GRBs 
may be observed by detecting the Cerenkov light emitted by high energy muons 
produced by neutrino interactions below a detector on the surface of the Earth 
(see for a recent review). The probability Pi/^ that a neutrino would produce 
a high energy muon in the detector is approximately given by the ratio of the high 
energy muon range to the neutrino mean free path. For the neutrinos produced 
in internal shocks, - lO^-* eV, P^^, ~ 1.3 x 10-*^(e^/lTeV) ||. Using (|^), 
the expected flux of neutrino induced muons is 

^P«10^^^km-V-^. (46) 

The rate is almost independent of e^b, due to the increase of P^^ with energy. 
The rate ( p6[) is comparable to the background expected due to atmospheric 
neutrinos p9| . However, neutrino bursts should be easily detected above the 
background, since the neutrinos would be correlated, both in time and angle, 
with the GRB 7-rays. A km^ neutrino detector should detect each year ^ 10 
neutrinos correlated with GRBs. Note, that at the high energies considered, 
knowledge of burst direction and time will allow to discriminate the neutrino 
signal from the background by looking not only for upward moving neutrino 
induced muons, but also by looking for down-going muons. 

The predicted fiux of ~ 10^^ eV neutrinos, produced by photo-meson in- 
teractions during the onset of fireball interaction with its surrounding medium, 
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Eqs. may be more difficult to detect. For the energy range of after- 

glow neutrinos, the probability P^^ that a neutrino would produce a high en- 
ergy muon with the currently required long path within the detector is P^^ ~ 
3 X 10-3(e^,/10i'^eV)i/2 This imphes, using (||), an expected detection 

rate of muon neutrinos of ~ 0.06/km^yr (over 27r sr), assuming fireballs explode 
in and expand into typical inter-stellar medium gas. If, on the other hand, most 
GRB progenitors are massive stars and fireballs expand into a pre-existing stel- 
lar wind, Eq. (^) implies a detection of several muon induced neutrinos per 
year in a Ikm^ detector. We note, that GRB neutrino detection rates may be 
significantly higher than derived based on the above simple arguments, because 
the knowledge of neutrino direction and arrival time may relax the requirement 
for long muon path within the detector. 

Air-showers could be used to detect ultra-high energy neutrinos. The neutrino 
acceptance of the planned Auger detector, ~ IC^km^sr Q, seems too low. 
The effective area of proposed space detectors |6^,0 may exceed ~ lO^km^ at 
eu > 2 X 10^^ eV, detecting several tens of GRB correlated events per year, 
provided that the neutrino flux extends to > 2 x 10^^ eV. Since, however, 
the GRB neutrino flux is not expected to extend well above ^ 10^^ eV, and 
since the acceptance of space detectors decrease rapidly below ^ 10"'^^ eV, the 
detection rate of space detectors would depend sensitively on their low energy 
threshold. 

Detection of high energy neutrinos will test the shock acceleration mechanism 
and the suggestion that GRBs are the sources of ultra-high energy protons, since 
> 10^'' eV (> 10^* eV) neutrino production requires protons of energy > 10^^ eV 
(> 10^^ eV). The dependence of ~ 10^^ eV neutrino flux on fireball environment 
imply that the detection of high energy neutrinos will also provide constraints 
on GRB progenitors. 

Inelastic p-n collisions may produce ^ 10 GeV neutrinos with a fluence of 
^ 10~^cm~^ per burst, due to either p-n decoupling in a wind with high neutron 
fraction and high, > 400, Lorentz factor (28|^, or to neutron diffusion in a wind 
with, e.g., strong deviation from spherical symmetry j?^. The predicted number 
of events in a Ikm'^ neutrino telescope is ~ 10yr~^. Such events may be de- 
tectable in a suitably densely spaced detector. Detection of ~ 10 GeV neutrinos 
will constrain the fireball neutron fraction, and hence the GRB progenitor. 

Detection of neutrinos from GRBs could be used to test the simultaneity of 
neutrino and photon arrival to an accuracy of ~ 1 s ('^ 1 ms for short bursts), 
checking the assumption of special relativity that photons and neutrinos have 
the same limiting speed. These observations would also test the weak equiv- 
alence principle, according to which photons and neutrinos should suffer the 
same time delay as they pass through a gravitational potential. With 1 s accu- 
racy, a burst at 100 Mpc would reveal a fractional difference in limiting speed 
of 10""'^^, and a fractional difference in gravitational time delay of order 10~^ 
(considering the Galactic potential alone). Previous applications of these ideas 
to supernova 1987A (see Q for review), where simultaneity could be checked 
only to an accuracy of order several hours, yielded much weaker upper limits: 
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of order 10^^ and 10^^ for fractional differences in the limiting speed and time 
delay respectively. 

The model discussed above predicts the production of high energy muon 
and electron neutrinos. However, if the atmospheric neutrino anomaly has the 
explanation it is usually given, oscillation to i^t-'s with mass ~ 0.1 eV plp7|j3^ , 
then one should detect equal numbers of I'^'s and Vr^s. Up-going r's, rather 
than /i's, would be a distinctive signature of such oscillations. Since i^^-'s are not 
expected to be produced in the fireball, looking for r's would be an "appearance 
experiment." To allow flavor change, the difference in squared neutrino masses, 
Z\m^, should exceed a minimum value proportional to the ratio of source distance 
and neutrino energy A burst at 100 Mpc producing lO^^eV neutrinos can test 
for Am^ > 10~^^eV^, 5 orders of magnitude more sensitive than solar neutrinos. 
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